Discrete Painleve equations, Orthogonal Polynomials on the Unit Circle 
and TV-recurrences for averages over U(N) — Pur and Py r-functions 



P.J. Forrester and N.S. Witte* 

Department of Mathematics and Statistics t(and School of Physics), 

University of Melbourne, Victoria 3010, Australia 

email: p.forrester@ms.unimelb.edu.au; n.witte@ms.unimelb.edu.au 

In this work we show that the N x N Toeplitz determinants with the 
symbols z M exp(-| y/i(z + l/z)) and (1 + z) M (I + l/z) v exp(te) - known 
r-functions for the P m / and Py systems - are characterised by nonlinear 
recurrences for the reflection coefficients of the corresponding orthogonal 
polynomial system on the unit circle. It is shown that these recurrences 
are entirely equivalent to the discrete Painleve equations associated with 
the degenerations of the rational surfaces — > E^p (discrete Painleve 
II) and — » (discrete Painleve IV) respectively through the 

algebraic methodology based upon of the afhne Weyl group symmetry 
of the Painleve system, originally due to Okamoto. In addition it is 
shown that the difference equations derived by methods based upon the 
Toeplitz lattice and Virasoro constraints, when reduced in order by ex- 
act summation, are equivalent to our recurrences. Expressions in terms 
of generalised hypergeometric functions q-^i n-^i are gi ven f° r the 
reflection coefficients respectively. 

f. Introduction 

There are now at least three approaches to systematically obtain recurrences for ran- 
dom matrix averages corresponding to r-functions for Painleve systems. One, exploited 
by the present authors is to use the theory of Schlesinger transformations within the 
r-function theory of Painleve systems. Another, due to Borodin ,;5], is based on a discrete 
analogue of the Riemann-Hilbert problem 0,13. A third due to Adler and van Moerbeke 
U is based on the theory of the integrable Toeplitz lattice and Virasoro constraints. In 
the work of the present authors, and of Borodin, there is an explicit connection with the 
discrete Painleve equations. Thus in all cases the recurrences for the r-functions involve 
auxiliary quantities which satisfy discrete Painleve equations. However the recurrences 
obtained in the work of Adler and van Moerbeke were not, in general, related to discrete 
Painleve equations. This then immediately raises the question as to the relationship be- 
tween the recurrences obtained by Adler and van Moerbeke, and the discrete Painleve 
recurrences. In this work, for the recurrences relating to r-functions for the Painleve III' 
and the Painleve V systems, we will answer this question by showing that in fact the 
recurrences obtained in |1] are transformed versions of the discrete Painleve equations. 
Moreover, we will show that the recurrences of Adler and van Moerbeke, obtained from 
their theory of the integrable Toeplitz lattice, follow from an approach based on the theory 
of orthogonal polynomials on the unit circle with semi-classical weights. In the situation 
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of a general weight the recurrence relations for the various coefficients appearing in the or- 
thogonal polynomial system are known as Freud or Laguerre- Freud equations |12l H]. This 
theory then provides a fourth approach to systematically obtain recurrences for random 
matrix averages corresponding to r-functions for Painleve systems. 

In this work attention will be focused on an orthogonal polynomial approach to r- 
functions for the Painleve III' and Painleve V systems defined as averages over the eigen- 
value probability density function for the unitary group U (N) with Haar (uniform) mea- 
sure (see e.g. |H| Chapter 2]), 

(1-1) (2 JnnI II l^-^'l 2 ' (z j :=e ie i,-ir<e j <n,j = l,...,N) 

^ ' ' l<j<k<N 

These r-functions are 

(1.2) r m '[N](t^):=t- N ^(l[ Z i 



N 

' U(N) 



(1.3) r v [N](t; /x, v) := < + *,)"(! + l/ztfe* 

\ A ■*■ / U(N) 

1 = 1 K ' 

where in 111.211 it is assumed (j£Z while in 11.31 it is assumed fi, v £ Z>o. By noting the 
identity 

(1.4) (1 + zY{l + l/zf = z (fl - v)/2 \l + zf +v 
we can rewrite 11.31 to read 

N 

(1.5) r v [7v]( i;M ,.) = (n^ ! " /2 ii+^rv^; 

\ A / U(N) 

which is well defined for 5t(/i + v) > — 1. Associated with 11.21 and 11.51 are the weight 
functions 

(1.6) z „JVt(,+^ 1 )^ z ( f ,-,)/2^ +z ^+u e t^ 

on the unit circle z £ T which have the special property of being semi-classical. On this 
latter point, analogous to the use of the term classical weight function for orthogonal 
polynomials on the line (see e.g. 0) we will call a weight function w(z) on the unit circle 
classical if its logarithmic derivative is of the form g(z)/f(z) with g(z) a polynomial of 
degree < 1, and f(z) is a polynomial of degree < 2. This gives w(z) = z^~"" 2 \l + z\^ +u 
as the only classical weight on the unit circle. The weights closest to classical with respect 
to the degree of the corresponding polynomials g and / are 

(1.7) z"e^ V ' t(z+z ~ 1 \ z ( ^- /)/2 \l + z\" + "e tz , 

and are to be termed semi-classical. In the case of the first weight in 11.71 with /i = 0, it is 
known that the corresponding orthogonal polynomials satisfy special recurrence relations 
[I3| [Tg] which lead to a recurrence for the corresponding U (N) average. We will show that 
this is also true of the general form of the first weight in 11.71 , as well as the second weight 
in (| 1 . 7^1 . thus giving recurrences for r ln [N] and r v [N]. In general these are different from 
those obtained in the Painleve systems approach, but rather coincide with recurrences 
obtained by Adler and van Moerbeke from their theory of the Toeplitz lattice and its 
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Virasoro algebra 4 . As already remarked, we are able to show that after appropriate 
transformations, the recurrences of the two approaches coincide. 

Another theme we wish to develop is the solution of the recurrences associated with 
11.211 and 11.51 in terms of generalised hypergeometric functions based on Schur polyno- 
mials. To define the latter let « = (ki, K2, • • • , kjv) denote a partition so that the 
non-negative integers with ki > K2 > ■ ■ ■ > kjv > 0, let s K (ti, . . . , tjv) denote the Schur 
symmetric polynomial, and define the generalised Pochhammer symbol 

N 

(1.8) [o]W :=H(a-j + l) Kj , (a)i:=a(a + l)...(a + l-l). 

3=1 

Also, with G k referring to a node in the Young diagram of « and a(i,j),l(i,j) the 

corresponding arm and leg lengths respectively (see [15] . define the hook length 

(1.9) h K = [a(i,j)+l(i,j) + l]. 

(i,j)GK 

With this notation the generalised hypergeometric series of interest is defined through a 
series representation [201 114] 



n in\ rt 1 )/ !. t i . \ I !]^ 1 ' ' ' s«(ti, . . . , ijv) 
(1.10) v F' q >{ai,...,a p ;bi,...,b q ;tx,...,t N ) = 2^171m TTTa) I 

for p,q € Z>o. The superscript (1) in p Fq indicates 11.1011 is a special case of a hyper- 
geometric series based on Jack polynomials and depends on a complex parameter d, the 
case d = 1 corresponding to 11.101 . 

In Section 2 we present some general formulas from the Szego theory of orthogonal 
polynomials on the unit circle [18] required for our study of Q 1 .2t and fEU i the later being 
carried out in Sections 3 and 4 respectively. In studying il l ,2t and l]1.5[l we first present 
the recurrence scheme following from the work of Adler and van Moerbeke and then 
proceed to show how the same schemes, and ones of lower order, can be derived from 
the theory of orthogonal polynomials on the unit circle. We make note of the solution 
of these recurrences in terms of the generalised hypergeometric functions, before revising 
the recurrence schemes following from the Painleve systems approach. The latter do 
not coincide with the recurrences resulting from the orthogonal polynomial approach, 
however we show that after an appropriate transformation of variables the recurrences 
are in fact equivalent. We label the particular discrete equations that arise according 
to the unambiguous algebraic-geometric classification of Sakai [17] by association with a 
degeneration of a particular rational surface into another, rather than the previous names 
employed. 

2. Orthogonal Polynomials on the Unit Circle 

We will consider orthogonal polynomials with respect to a complex weight function 
w(z), analytic in the cut complex z-plane. The latter specification means w(z) possesses 
a Fourier expansion 

(2.1) w(z) = ^2 w k z k , w k = / ~^j^ w ( z ) z ~ k > 
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where T denotes the unit circle \z\ = 1, appropriately deformed so not to cross the cut, 
and z — e lS , 9 6 (—71", 7r]. For e = 0, ±1 we define the Toeplitz determinants 



(2.2) r N [w] := det 



dz , N e-i+fc 



-w(2)z e 



= det [io_ e+j _ fc ] < J . <fc < JV _ 1 . 

0<j<fe<JV-l 



27T12 

In the case e = 0, and T the unit circle without deformation, by virtue of the identity 



(2.3) det 



dz , s e+j-fc 



2niz 



w(z)z e 



N 



TT zfw(zi) 

0<j<k<N-l ' t=l ' " M 



we see that 11.211 and (11.51 1 . in the cases (igZ and (J,,v £ Z>o at least, can be expressed 
as /jv[w] with w(z) as in 11.71 . In certain circumstances the weight is real and positive 
w(z) = w(z) and thus the Toeplitz matrix is Hermitian, w/t = w-k, but in general this 
will not be true. 

If -Zjv[w] is non-zero for each N = 1,2,... then there exists a system of orthogonal 
polynomials {<f> n (z),n = 0, 1, . . .} with the orthonormality property 

f dz 

(2.4) J ■^—r-w{z)4> m {z)(l) n (z) = <5 m ,„ 

where 

(2.5) w(z) := HM. 

Wo 

We introduce special notation for the various coefficients in (j) n {z) according to 

n 

(2.6) <frn(z) — K n z n + InZ 71 ' 1 + m n z n ~~ 2 + . . . + <f>n(0) = y~] C n jZ j , 

3=0 

where without loss of generality k„ is chosen to be real and positive. We also define the 
reciprocal polynomial by 

n 

(2.7) z n $(l/z) ^cnjz"-', 

3=0 

where c denotes the complex conjugate. A fundamental quantity is the ratio r n = 
4> n (0)/tt n , known as reflection coefficients because of their role in the scattering theory 
formulation of orthogonal polynomial systems on the unit circle. 

From the Szego theory |18| these coefficients and their complex conjugates are related 
to the above Toeplitz determinants 1)2. 3|l by 

(2.8) rjv = ( _iyvfjH fjv = ( _ 1)W ^H 



I%H' " J I°n 



w 



In the case that w(z) is not real, fjv (notwithstanding the notation) is not the complex 
conjugate of rjv but rather an independent variable. Note that the normalisation 1)2.5^ 
implies that ko = 1 and thus ro — fo = 1. Knowledge of {vn}n=o,i,..., {fiv}.iv=o,i,... is 
sufficient to compute {/jv[iu]}jv=o,i,... • For this one uses the g enera l formula |18| 

(2-9) {I° N [w}) 2 =!-™- 

Our fundamental task is then to obtain recurrences determining the tm and fjv. 
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For this purpose we require further formulae from the Szego theory. First, as a con- 
sequence of the orthogonality condition we have the mixed linear recurrence relations for 

(fin and 0*, 

(2.10) K,„(j) n+ i(z) = Kn + xZ<fi n {z) + (fi n+ i(0)(j) n (z) 

(2.11) (z) + 4>n+l(0)z<fin(z), 
as well as the three-term recurrences 

(2.12) K n (j} n (0)(j> n + 1 (z) + K n _l0 n + l(O)z0 n _l(z) = (k„0„ + i(O) + K n + 1 (j>n(0)z)(j} n (z) 

(2.13) K n 4> n {Q)(f)* n+1 (z) + K n -l(f>n + l(0)z(j} n _ 1 (z) = (K n (j) n + l(0)z + K n + 1 (j) n (0))(j)n{z). 

From the latter one can derive the analogue of the Christoffel-Darboux summation formula 



(2.14) X>j(*)&(0 



K( Z )<t>n(0 ~ <<Mz)<MC) 



1 - * 



/cy -. r\ _ ^n+lfcX+lCC) ~ <Pn + l{z)<t>n + l(Q 

1 - zC 

for zC, 7^ 1. Identities from the Szego theory that relate the leading coefficients back to 
the reflection coefficients are 

(2.16) 4 = K Li + l<Mo)| 2 , 

/ 1 

(2.17) ^ = Y J v j ^f h 

Kn 3=0 

n — 1 , 

(2.18) = 2^ r j+irj -! + r j+irj 



3=0 



Another relevant formula is 



N-l 



(2.i9) i° N [w] = n Kj 

3=0 

Finally, with n n denoting an arbitrary polynomial in the linear space of polynomials with 
degree at most n, we can check from the structure 12. 61 that 



z<fi„(z) = 



Z (fin(z) 



In 


ln + 1 \ , 




Kn + 1 J 


In 


\ ! m n 


Kn / 




In 


ln+2 



Kn+1 

In ( + 1 In \ TTln rZn 



<fin~l{z) + TT„-2 



(2 .20) + \ h±L ( ±±1 - k.) + r^L_!^±iU n(0 ) +7Tn _ 1 

[ Kn + 1 \rtn + 2 K„ J K,„ K n + 2 J 

<f>'n(z) = U K " (fin-l(z) + 7T n _ 2 
Kn—1 

Zlj>' n (z) = n<f>n(z) ^— (fin-l(z) + 7T„_2 

fan— 1 

z 2 <fi' n (z) = w K " n+ i(z) + I (n — 1)— — n 1 </!>„(z) + 7r n _i 
where ' denotes the derivative with respect to z. 
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How we use the above formulae to produce recurrences for rjv and fjv differs for the 
orthogonal polynomial systems corresponding to the two different weights in l|1.7[l . Let us 
then treat the two weights separately. Let us refer to the order of a difference equation in 
r„, f n as q/p where q 6 Z>o refers to the order of r n while p G Z>o refers to the order of 

fn- 

3. The P m / System 
For non-integer values of p the weight function 

(3.1) w(z) = z li ei y/i( - z+ '~ 1) , 

has a branch point at z = 0. Cutting the complex plane along the negative real axis 
(—oo,0] we deform the contour in 12. II to the contour starting at — oo, running along the 
real axis on the negative imaginary side to z = — 1, following the circle \z\ — 1 in the 
anticlockwise direction to return to z = — 1 on the positive imaginary side, then returning 
to — oo along this side of the negative real axis. This contour is standard in the theory of 
the Bessel function (see |19l.pg. 363). Denoting this contour by C, and noting that the 
integral representation of the Bessel function of pure imaginary argument gives 

(3.2) / e £»M='-<^. 
we see that 

(3.3) In[w] = det[7 M+e+ j_ fe (v^)]j,fc=i,...,jv- 
For general fi we then define 

(3-4) T m '[N](t;iJ.) = det[I ll+j - h (y/i)] j , ll=1) ..., N , 

which is consistent with 11.21 in the case fi G Z. 

Recurrences for the reflection coefficients rjv, f jv in the case of the weight 13.11 can be 
deduced from the work of Adler and van Moerbeke flj. In their case 3, one specialises their 
parameters d\ = = f[ = 72 = 7" = 72' = 0, 7 = fj, and one sets Pi (2) = Pi[z) = 1 /2\ftz. 
Then making the identification rjv = xn, tn = J/jv it follows from jl], eq. (0.0.17) that 

(3.5) 1 /2Vtv N {r N +i +rjv_x) +Nr N = 0, 

(3.6) V 2 Vtv N (f N+1 + fff-i) + Nfff = 0, 
where vn ■= 1 — rjvfjv. After specifying the initial conditions 

(3.7) r =r = l, ri = — r , >, , ^ = — A r , ' , 

which follow from 12.81 and 13.31 . the recurrences 13.51 and 13.61 uniquely determine 
TNi 'Tn f° r N = 2,3, . . .. We note that the order of 13.51 and 13.61 is 2/0 and 0/2 re- 
spectively and the parameter p, does not appear explicitly. In addition we observe that 
13.51 and 13.61 have the familiar form of the discrete Painleve equation associated with 
degeneration of the rational surfaces Dg 1 ' — > |17| (discrete Painleve II). We now seek 
a derivation of 13.51 and 13.61 using the formulae of Section 2 specialised to the weight 
l|3.1|l . In addition we will show that the orthogonal polynomial theory can be used to 
derive a pair of coupled difference equations for rjv,fjv, both of order 1/1. It will then be 
shown how I|3.5I3.6|I can be deduced from these equations. 



Proposition 3.1. The reflection coefficients f2.fip corresponding to the Toeplitz determi- 
nants satisfy the coupled 1/1 order recurrences 

(3.8) 1/2 Vi(r N+1 f N + rjvfjv-i) + N - rNrN _ /i = 0, 

1 - rjvrjv 

(3.9) V 2 Vt{f N +ir N + r N r N -i) + N - rNrN _ +/i = 0, 

1 — r N r N 

with the initial conditions \3. 7| ). 

Proof. We adapt a method due to Freud |12|. and consider two different ways to evaluate 
the integral 

(3.10) Ji := / -^-z 2 w' (z)4> N (z)<j> N+1 (z). 

Integrating this by parts, employing I2.20H and the orthogonality conditions shows 

(3.11) Ji = — (AT + 1)— ^ — \-(N + 1)— ^ii. 

Kn+1 Kjv 

Alternatively we note from 13. Ill that 

(3.12) ^ = £ + x/ a Vt(l-4), 

w z z 

substituting this in 13.10L one can again use I2.2UH and the orthogonality conditions to 
show 

(3.13) Jl=M ^L + i/ 2 ^f-^- lN+2 Kn ). 

KN+l \KN+1 KJV+2 KJV+i J 

Equating 13.111 and 13. 1311 and eliminating In using (12.1711 gives 

EHJ- T° deduce lEHI 

we apply an analogous strategy to 

(3.14) J 2 := J ^iu'(*)^+i(*)M«). 

□ 

To deduce <EHJ from ^3, J3~^ we first subtract (££HJ| from O to find 

(3.15) V 2 Vt(LN+i + L N ) + 2/Lt = 0, Ljv ~ fjvrjv-i - rjvfjv-i. 

Using 13. 71 and a Bessel function identity this equation is to be solved subject to the 
initial condition 

(3.16, I,—* 
It follows that the solution of 113. 1511 is the constant 

(3.17) L N = -^, N = 1,2,3,... 

vt 

and thus 

(3.18) V 2 Vt(f N +ir N - rjv+ifjv) + ^ = 0. 

Using this to substitute for [i in Ij3,8^ and l|3.9[l gives i'.'>.~>t and l|3.9[l respectively. Fur- 
thermore l|3.18fl can be summed to give 

(3.19) i/ ay /i(hL-hL) +ll N = o. 

\k n kn J 



s 



We turn our attention now to formulae for t hi [N], rjv and rjv in terms of generalised 
hypergeometric functions. From earlier work |14) . we know 

(3.20) T UI '[N](t;ri = (^J II r(/+ p) (i N + w * 1 ' - ■ - ' ^)[ tl= ... = t JV= t/4. 
It follows from this and f3~3l. and <2~%l that 

,-. v = (-1)"/^* 1 o^i (1) (;^ + l + M;ti,...,tiv)| 



2/ (m + 1)jv F x (1) (;iV + /Lt;ti,...,t N ) lti=...=tjv=t/4' 



(3.22) Fjv = (-1)" ( I ( M ) JV , 

Note that the small-t expansions are more forthcoming from these formulae than from the 
Toeplitz determinants. 

Our final point in relation to r [N] concerns the known recurrence scheme for r m [N] 
in terms of the variables pjv, qjv specifying the corresponding Hamiltonian in the Painleve 
systems approach to Pup . The idea here is that we start from the Hamiltonian for the 
Pup system 

(3.23) tH™ = qlpi - (qi + viq n - t)p n + l /i (vi + v 2 )q n , 

with special parameters (1)1,1)2) = (v[°\v^) = (p,—p,) and corresponding special values 
p — po = 0, q = qo for some particular go (see e.g. \W\)- A sequence of Hamiltonians is 
constructed from this seed by application of the Schlesinger transformation Ti with the 
action on the parameters Ti • (1)1,1)2) = (t>i + 1, V2 + 1) and some explicit actions on p and 
q involving rational functions of p\( vi ,v 2 ) an d Q\(v 1 ,v 2 )- Thus we set 



(3.24) tH™' := tH U1 ' 



rr(4 0) .4 0) )' 



and introduce the corresponding r-function t™ by the requirement that 

(3.25) H™' :=±lo gT ™'. 

We know from that with (v^\v^) — (p, —p) the sequence {t™ }n=o,i,2,... is realised 
by 

(3.26) t™' =t- n » / \ m '[n\{t ]l A . 

I ti— >4t 

With p n := T"p\ (o) (o). and q n := T"q\ ( ) <o), further theory associated with Ti led 
to the following recurrence scheme. 

Proposition 3.2 (0). We have 

T^[JV + l]T m '[JV-l]| _ r)vr _ 10 x 



(3.27) 



( T ni'[jv])a 



lt>->4t 




PJV+1 = 


T (piV 


9iv+i = 


t 



CS.2M Pjv+ i = ^f(pjv-l)-M^+l (iV = 0,1,...) 

<•■'>■*» vm,--— + f / JV+ ^ 1,, (iV = 0,l,...) 

gjv gjv [gjv (pjv — 1) — + * 
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subject to the initial conditions 



(3.30) Po = 0, q =t^-lo g r^ 2 I^Vt)\ 

at lt>->4t 

(3.31) r m '[0] = l, r lu '[l]\ t ^ it =I,(2Vt). 

Our interest is in the relationship between the variables pn,<1n in these recurrences 
and the reflection coefficients of Prop 13. ll This is given by the following result. 

Proposition 3.3. The Hamiltonian variables qM,PN in Prov \3.°A are related to the reflec- 
tion coefficients in Prov \3.f\ by 

(3.32) pjv = 1 — r N f N \ t ^4t, 

(3.33) 9JV = _VI^| . 

TjV 1*1— >4t 

Proof. The equation 13.3211 follows immediately upon comparing 13.271 1 with 12.91 . Sub- 
stituting 13.321 in 13.281 we see that 13.181 results if we also substitute for according 
to 13.331 . Furthermore, if we combine 13.281 and 13.291 into 

t N 

(3.34) q N + = — , 

qN-i Pn 

we see that the substitutions 13.331 and 13.321 in this equation yields 13.51 . □ 

Thus the structure of the recurrences is such that the transformation equations l|3.33^ 
and l|3.32fl can essentially be determined by inspection. However when we come to study 
the analogous recurrences for the Pv r-function 11.51 this is not possible and a more 
systematic procedure is called for. With this in mind, let us then present a more systematic 
approach to the derivation of l|3.33fl . 

For this purpose, in addition to the shift operator Ti which increments N in r m [N], 
we introduce the other fundamental Schlesinger transformation of the P n i' system T2 with 
the action on the parameters T2 • (vi,V2) = (vi + l,Wa — 1). Recalling that in t™ [N] 
(v^, «2°') — (Mi — A 4 ) i this operator then increments fj, by unity. From |1(J| we know that 
for general parameters 

(3.35) Ti-tH 111 ' =tH m ' -q(p-l), T 2 ■ tH m ' = tH™' - qp, 
and thus in particular 

(3-36) q N = t— log p TN 

at 1 2 ■ tn 

(3.37) =1/2 (AT- fi) -t-| log K 2 N r N \ 

at Ui-»« 

where to obtain the second equality use has been made of 13.251 . 12.81 and 12.191 . To deter- 
mine the ^derivatives of the orthogonal polynomial coefficients we find the ^derivatives of 
the polynomials themselves. Differentiating the orthonormality condition 12.41 with w(z) 
given by 13.11 we have 

(3.38) = -— S mn + [ -^—w[j) m J r 1 /2Z(t> m \4>n+ [ -f^w4>m\^ n + X h Z(j) J, 

Wo J j 2mz J T 2niz 
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where ' represents differentiation with respect to t. Now <f> n + 1 /2 z<f> n is of degree n + 1 
in 2 and has no components in 4> m for m > n + 2. With this established it follows from 
113.3811 that 4> n + V 2 z 4>n has no components in cj> m for m < n — 2 either and so 

(3.39) cj>„ + 1 /'2 z<j) n = a n <j}„ + i + b n (j) n+ i + Cn4> n -i. 

Equating coefficients of the highest monomial in 13.391 and recalling 12.61 gives a n — 
fv n /2K„+i while 13.381 in the case m = n — 1 yields c n — —n n -i/2K„. Finally, setting 
m — n in l|3.38|l and recalling 1)3. 2|l shows 

In 



(3.40) 



b n + b n = 



Substituting for a n and c n in l|3.39|l . and using the three-term recurrence l|2.12^ to substi- 
tute for a n (j>n+i 



(3.41) 



1 /2 z(j>„ we deduce 
r„+i 



1 H z 



Recalling 12.61 it follows that 



(3.42) 
(3.43) 



— y- + V 2 (r n +ir n + r„+ir„) 



V 2 ( r n^ 



r n -i 



(3.44) — = V 2 (r n +i - r n -!) \ — 

Making use of the first two of these relations in 13.371 and simplifying reclaims 13.331 . 

4. The P v System 

The weight function 

(4.1) W(2) = (1 + 2)"(1 + 1/V)V Z , 

is analytic in the cut plane z G C\(— oo,— 1]. For 5R(/i + v + 1) > the singularity is 
integrable at z = —1 so with this restriction there is no need to deform T in JO)- The 
transition from the Pv to the Pur system can be achieved by making the replacements 



(4.2) 



v - ii, 



2vzi 



t/Av 



with t, N, n fixed and then taking the limit v — > oo. In this way the Toeplitz determinant 
11.31 reduces to 11.21 . In |U] we have shown that with w given by 14.11 



(4.3) 



dz 
2niz 



n w{z) 



r(/x + 1 - n)V(u + 1 + ; 



iFi(—n — v; fx 



l;-t), ne. 



and thus ljv[w] as specified by 12. 2t can be made explicit. In particular it follows from 
tOl . ll2~2l and JTH) that 



(4.4) 



ri 



iFi(-i/+l;/i + 2;-t) 



A* 1F1 



l;fi;-t) 



As with the weight 13.11 . recurrences for the reflection coefficients rjv,fjv in the case 
of the weight 14.11 can be deduced from the work of Adler and van Moerbeke Thus 
in their first case one specialises their parameters 7" =72 =72 = ^2 = 0, d\ = —1,7 = 
— f,7i = fi + v which gives their (a,b,c) = (1, 1,0), and one sets -Pi (2) = tz,P?.{z) = 0. 
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With the identification xn = rjv,yjv = rjv, from [I] eq.(0.0.14) we then read off the 



recurrence 



(4.5) — trN+irN + irjvrjv-i + (N + 1 + v)tntn+i — (N — 1 + v)rN-i?N = 0, 
which is of order 2/2. Also [1] eq.(0.0.15) reads 

(4.6) v N +i [N + 1 + v + tr N+2 f N ] - v N [N + v + tr N+ ifN-i] + r N +ifN 

= ni(l + v + tr'/fo) + rif , 

where vn '■= 1 — rjvfjv which is also of order 2/2. The right hand side of l|4.6^ can be 
simplified. First we note that it can be written in terms of the Fourier components of 
14.1|l according to 



(4.7) (y + 1)(1 — w-iwi) + t(ni-i — w-2) — w-i 

and thus in the light of il l.MI in terms of the confluent hypergeometric function. Utilising 
the contiguous relation for the confluent hypergeometric function i_Fi(a; b; —t), [5j (13.4.7, 
p. 507) in the two cases a — —v, b = fj, + 1 and a — — +1,6 = fi + 2 one can show this is 
precisely unity, and so l|4.6^ simplifies to 

(4.8) ujv+i [N + 1 + v + tr N+ 2f N ] - v N [N + v + trjv+ifjv-i] + rjv+ifjv = 1. 

We can use an orthogonal polynomial approach to derive recurrences for rjv , f jv , and 
from equations used in the derivation f4.50 and l|4.8|l can be reclaimed. 

Proposition 4.1. The reflection coefficients rN,r~N for the orthogonal polynomial system 
on the unit circle with the weight U-lf satisfy the coupled system of 2/1 and 1/2 order 
recurrence relations 

(4.9) (JV — 1 + v)r N -i + {N + fi + t)r N = -t(l - r N f N )r N+1 + tr%f N -i, 

(4.10) (N + 1 + u)f N+ i + (TV + fi + t)f N = -t(l - r N f N )f N -i + tf 2 N r N+1 , 
subject to the initial conditions \4-4f - 

Proof. As with the derivation of the recurrences of Proposition lM.il we will take the Freud 
approach, but now working with two integrals for each recurrence rather than the single 
integral for each recurrence required to derive the recurrences of Proposition 13.11 The first 
integral we consider is 

(4.11) J i: = / ^-{1 + z)w\z)4>n{z)J^{z). 

Integrating this by parts, employing I2.20H and the orthogonality condition shows 

(4.12) Ji:=N hL-^ N + i) l Ji+L, 

KJV ACJV + 1 

On the other hand we note from 1)4. l|l that 

(4.13) — = T^--+^ 

w 1 + z z 

substituting this in 14. lit , an analogous calculation shows 

(4.14) Jl = il - v (hL-h±L\ +t ( l JL-k+L\ +t 
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Upon employing 12.17|l one can equate Q4.120 and l|4.14[l and solve for ijv+i to obtain 

(4.15) 1n+1 = -fi - t(l - f N r N +i) - (A + v)r N f N +i. 

Note that one could perform another differencing at this point and eliminate In in favour of 
TM,rN but then /i would disappear from the ensuing relations. This is a clear indication 
that the recurrence system would be raised unnecessarily in order, so we seek another 
relation for In, and this is found by considering the following integral 

(4.16) J 2 := J — z(l + z)w'(z)<j> N (z)<t> N+1 (z). 
The methods of evaluation used to derive 14.121 and 14.141 now yield 

(4.17) j 2 = (i v + l)W _ (Ar + l) J^_^v ± i ) 

and 

1oN , kn , . ( km In In+2 k-n 

(4.18) J 2 = fx h t h 



KN+l \Kn+1 Kn+1 KN+2 Kn+1 , 

respectively. Equating 114. 1711 and 14.18L and again employing 12.17L it follows 
(4.19) =N+l-[N + l + ii + t- t(r N+2 f N +i + r N +if N )] (1 - r N+1 f N +i)- 

KJV+1 

Equating 14.1511 and 14.191 1 gives 14.911 . The second recurrence can be found in a similar 
manner by eliminating In+i from expressions arising from evaluation of the integrals 

dz 



■h ■= I -^r^(l + z)w'(z)4in+i(z)<Pn(z), 



and 

dz 



Ji~ I - . z(l + z)w (z)4>n(z)4>n(z). 
Jrj. 2mz 

The two expressions for Zjv+i are respectively 

(4.20) }n+i_ _j_ i) rjv+1 r JV+1 _ [y _|_ tr N+2 r N )(l - r N +if N +i) 

= -v-(N + fi + t)r N +ir N - t(l - r N f N )rN+irN-i 

(4.21) - t(l - r N +ifN+i)r N +2fN + tr 2 N+1 f 2 N , 

and this yields 14.101 . □ 

Let us now show how we can derive the recurrences of Adler and van Moerbekc 14.51 , 
14. 81 from the workings of the proof of Proposition 14.11 As already remarked, subtracting 
from 114. 1511 the same equation with A replaced by A — 1 and recalling 12.171 gives a 
difference equation in rjv, fjv. This difference equation is in fact precisely jHH - Note that 
with Ajv the forward difference operator with respect to A, 1)4.15(1 itself can be written 

(4.22) A N [{N + v)— -t— ] = -/!- 1. 
Summing this over A shows 

(4.23) (N + u)!Z--tlZ-=-(fi + t)N, 

which then is the basic relation underlying l|4.5|l . The recurrence l|4.8fl is obtained by 
subtracting from l|4.20|l the same equation with A replaced by A — 1, and recalling |2.17p . 
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We consider next formulae for t [N](t; p,, v), rjv, rjv in terms of generalised hyperge- 
ometric functions. From earlier work |14| . |11) we know 
(4.24) 



t-i = . . . — t !\r — — t 



Noting from FFZl . llOl and itTll that for the weight lETTt 

(4.25) I 1 N [w] = T V [N](t;v + i^M = T V [N](t; + 
it follows from this and 12.81 that 

(4.26) r„ = (-!)", ^ 1 F 1 W(-,+ l ; JV + l + W t 1 ,-,t Jf )| 

(M + 1)jv ^^(-i/jJV + ^ti,.-.,^) Iti=...=tjv=-t 

(4.27) fjv = (-l) r 



, A (p) N ^^(-v-liN -l + n;ti,...,t N )\ 



tl=. ..—t AT — — t 



(^ + 1 )jv 1 F-t 1) (-v;N + f i;t 1 ,...,t N ) 
An immediate consequence of I4.26H and 14.27H is the explicit values at t = 0, 

(4.28) rN = { - 1 f(^, r N = (-lf Mn 



(4.29) 



(m + 1)jv' v y (fH-l)w' 

l N __ Nv In Np 
kjv N + fi' km N + v 



Again we see that the transition from the generalised hypergeometric function in the 
Pv system 14.241 to that of the P n i' system 13.201 is facilitated by making the replace- 
ments I4.2|l and taking the limit v — > oo 



(4.30) x FP{a-v\N + \i;tx,...,t N ) 



t\=...=ipj=—t/Av 



F^ ) (;N + p;t 1 ,...,t N ) 



\t!=...=t N =t/4 

for all fixed a, N + p, t. Thus the corresponding limiting forms for the reflection coefficients 



V i iV \ III' -V | VM -III' 



are 

(4.31) r\ ^ (^) r\". r\ - 1^) A" . 

where we distinguish the two systems only when some confusion could arise. Another 
consequence of this limiting process is that the recurrences in r jv , f n lETTrJl reduce 

to a sum of 13.81 and (13. 1811 . and of 13. 9> and 113. 1811 respectively. In addition we find that 
l|4 r 2ljl reduces to fTHty . 

It remains to compare the recurrences Ij4.9^ . I|4.1()^ determining t [N] through the 
recurrence J2jj, and recurrences satisfied by the Hamiltonian variables gjv,pjv in the 
Painleve systems approach to Pv which also determine i~ v [iV] through a recurrence 
In relation to the latter, let v\ + v 2 + V3 + V4 — and introduce the Hamiltonian 

(4.32) tH v = q(q — l)p(p + t) — (v 2 — V\ + v 3 — Vi)qp + (v 2 — Vi)p + (vi — V3)tq. 

Starting with a special solution when V3 — V4 — V3 — = 0, q = qo = 1, p = po for 
some particular po (see 0), a sequence of Hamiltonians is constructed by application of the 
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Schlesinger transformation T 1 with the action on the parameters T 1 • (vi, v?, V3, 1)4) 
(vi — 4 /4 , V2 — V 4 > v ?> ~ V 4 > v i + 3 / 4 ) to obtain 



(4.33) 



l^ji — >(vi — n/ 4,V2— Ti/4,TJ3 — n/4,^4+3n/4) 



The corresponding sequence of r-functions t„ are specified so that 
(4.34) 



rrV O , V 

H „ = — log T„ . 



We know from [5] that with v}' 
{TvT}n=o,i,... is realised by 



(o) „,(<>) 



dt 

(0) (0) 



(0) (0) 



the sequence 



(4.35) 



rv + ^ + i) 



J] I> + Z + l)r V [n](t; M ,*,). 



Furthermore {r v [iV]}jv=2,3,... is determined by the following recurrence scheme [U]. 
Proposition 4.2. Let 

(4.36) Xjv = {pn + t)q N + 1/2 ^i, j/jv 
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The sequences {T V [-/V]}jv=o,i,..., {xn}n=o,i,...> {yjv}jv=o,i,... satisfy the coupled recurrences 



r [N + l]r v [N- 1] 



(rV[AT])2 



IJV + CEjV-1 = 



(4.37) (N + v)- 
(4.38) 

(4.39) Vnun+i — t 

subject to the initial conditions 



Xn v — V2 fj, 

ytt 

t N 



■N 



yjv 1-vn 

X N + N + l + U+ 1 /2fl 



(4.40) 
(4.41) 



x =t + 1 /2 n + t— log 1 F 1 (-v: n + 1; -t), y = 1 



By adopting a strategy analogous to the method of derivation of 113.3311 . 113.3211 given 
below the proof of Proposition 13.31 the relationship between the variables XN,yN and 
rN,r~N can be deduced. 

Proposition 4.3. The Hamiltonian variables xn , yjv and the reflection coefficients rjv, f\zv 
are related by the equations 



(4.42) 



xn v — V 2 A 4 

J/jv 



XJV 1/ ~ 1 /2 fl 

VN 



1 + 



+ iV= (JV + !/)(l-nvrjv) 



N 



(4.43) 
(4.44) 



t (1 - ^ivrjv) 



(1-Vn) 



a^jv h V 2 A 4 



+ JV = -t^ i (l-r w f w ) 
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Consequently 

(4.45) 

(4.46) 



Vn = - 



v + tr N+1 r N 



x N 



rN[tr N +i + (N + v)r N ] 

t , . 
1 /2 fi = —tr N+1 r N . 

VN 



Proof. The equation 14. 421 follows immediately upon substituting for r v [N + l]r v [7V — 
1]/(t [iV]) 2 according to the right-hand side of the general relation 12.91 in 14.371 . For the 
remaining two equations, adapting the method of derivation of 13.331 . 13.321 given below 
the proof of Proposition ^. 31 we require a shift operator that has the action /i i— > + v i— > 
v — 1 and leaves N fixed. Such an operator is the Schlesinger transformation r 2 _1 with the 
action on the parameters r 2 _1 • (vi,t>2, 133,1)4) = (vi + Vi — 1 /a , V3 — x /4 , V4 — 1 /a ) which 
was studied in |1U|. From the explicit form of T 2 _1 in terms of operators associated with 
the root lattice ^3, and the actions of these operators on the Hamiltonian and associated 
variables in |10| . we can compute that 



(4.47) 



tH v =tH v + [q 



Vl — V S 



Vi — 1)4 



(4.48) T 2 -tH v = tH v + ( g -l)( p -H^l 



Recalling KMl . E"351 and llOrJll it follows that 
(4.49) 



1 4- Vl - Vs 
p 



t - (V3 — V4) 



+ t + V 3 —V4 



XN V — V2 (1 

VN 



N + v 



\ XN ~ W ~ V 2 m) (! ~Vn)- ' 



-N = t-lo g r n 



Xn h V 2 A 4 

VN 



d 



+ N = t-lo g r n . 



(4.50) 0—Vif) 
To find the f-derivatives we differentiate the orthonormality relation 



(4.51) 



: 



w x 

Omn "T 



dz 



W[(j) m + Z(f>r, 



+ 



dz 



— W(f>m4>r, 



Wo J 2mz J Ziriz 

and the case n < m — 2 indicates that (f> m + z(j> m has no components in <j> n so 
(4.52) <f> m + z<j} + 

Consideration of the coefficients of the highest power in z gives a m = K m /K m +i, and the 
case n — m — 1 shows e m = 0. The t-derivative of the orthogonal polynomial is 



(4.53) 

and this implies 
(4.54) 



d n + 



ft n n+ i(O) 



Kn-l<ftn + l(0) 
K n + l<t>n{0) 



Z(j>n 



r„+i ( r, 



Next we examine the case n > m + 2 and find 

(4.55) (j> n = b n <j> n + c»^„_i. 
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Taking the case n = m + lwe infer c n = — K n _i/K n . The t-derivative of the orthogonal 
polynomial can then be written as 

(4.56) <f> n = b n <f) n — 1 4> n -i, 
from which we deduce 

(4.57) f JL = - fn _JL-, 

Substituting 14.541 . 14.571 in 14.491 . 14.501 gives 14.43I4.44H respectively. 

To derive 14.451 . multiply both sides of 14.431 by (1 — yjv)/2/jv, then substitute for 

y l ~VN v t 

X , X := Xn v — V 2 Mi 

VN UN 

using 14.421 to deduce 

[(TV + „)(i - r N f N ) — N] (l l —) N^M. = t r Jt±l {1 - nvr»i^. 

Solving this equation for yN gives 14.451 . For the derivation of 14. 461 we write 14.421 in 
the form 

(X + v){l - y N ) — —{N + v)r N f N y N + v. 
Substituting 14.451 for yN and simplifying we obtain 14.461 . □ 

The Hamiltonian variables qN,PN or xn,Un in the Py theory go over to those in the 
Pin/ theory under the replacements 14.21 and upon taking the limit v — > oo 

Ia col V III' V III' 

(4.58) q N -> 1 - p N , Pn -> Qn 



v — »oo v — >oo 



/A rn\ V , / HI' /l V 1 

(4.59) x N -> V 2 M + 9iv (1-Pjv )> Vn -> 7777- 

V — 'OO v — >oo — PjSf 

Using these transitions we find that the recurrence relations <Ofl . K3§i reduce to tSTM 
upon using lEHTl (which follows from lTH9l and KFM ). and to J-TSSt respectively. 

In addition to the formula 14.451 for yjv, we can obtain a different formula by making 
use of 14.441 . This allows a pair of 1/1 order difference equations for rjv, fjv to be deduced, 
thus reducing the 2/1 system of Proposition 14.11 down to the same order as the coupled 
system 14.381 . 14.391 satisfied by xn,Un- 

Theorem 4.1. The reflection coefficients satisfy the coupled 1/1 order recurrence relations 

(4.60) (1 - r N f N )[tr N+ i + (N + v)r N ][tf N -i + (N + v)f N ] 

= [{N + v)r N f N + n][v - (N + u)r N r N ] 

(4.61) t 2 r&f&_i +t(v-n- t)r N f N -i - (N + u)(N - 1 + v)f N r N -i 

— {N — 1 + v)tr N -if N ^i — (N + u)tr N fN — fiv = 

Proof. Multiplying both sides of 14.441 by l/yjv and substituting for X(l — yN)/yN using 
14.421 then solving for j/jv shows 

N + + v + i%^(l - r N f N ) 

(4-62) y N = + (M 1\ - ■ 

fj, + (N + v)r N rN 

Equating this with Ij4.45^ and solving for rjv+i gives l|4.6()^ . The second recurrence follows 
by eliminating rjv+i between IJ4.6O0 and l|4.9|l . □ 
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As with the other recurrences for rjy, fjv in the Pv system, we find l|4.6()|l . I|4.ril|l under 
the replacements l|4.2|l and taking the limit v — + oo assume the forms of the recurrences 
in the P m i system (13. 81 . 13.181 respectively. 

5. Concluding Remarks 

In our work in addition to the U(N) averages (11.21 . 11.31 . an TV-recurrence was 
also obtained for 

N / 1 \ M 

(5.1) r VI [iV]( i ; M ,a; 1 ,^; ? ) = (n(l-Cx[l ) _^))^ 2ei ! 1 + ^! 2 " 1 (^J (i+^H^. 

where Xj = 1 f° r S J, Xj — otherwise. As the notation suggests, this is a known 
r-function for a Pvi system [TTj. Both r HI [TV] and t v [TV] can be obtained as degenera- 
tions of (15. It . or equivalently the weights (11.711 can be obtained as limiting cases of the 
"master" semi-classical weight function underlying 15. II . For general parameters, a recur- 
rence scheme based on the discrete Painleve equation associated with the degeneration 
of the rational surface — > (discrete Painleve V), was given as a consequence of 
the Painleve system theory of Pyi (recurrence schemes for special cases of the parameters 
have also been given in 0, 0). It is also true that an TV-recurrence for 15.11 can be 
deduced from the Toeplitz lattice approach of Adler and van Moerbeke But here we 
have found that in the cases of the U(N) averages (|1.2|l . l|1.3|l the approach of 0] leads 
to equivalent results as do those obtained from an orthogonal polynomial approach. One 
therefore suspects the same will be true in relation to l|5.1fl . and that furthermore the 
corresponding recurrences are transformed versions of the discrete Painleve V equation 
found in This is indeed the case, but the details do not fit well with the scheme of 
the present paper (in particular the Freud approach to the recurrences for r„,f n is now 
inadequate) so will be reported elsewhere. 

Another point of interest relates to the TV-recurrences for Hermitian matrix (as op- 
posed to unitary matrix) averages in which the weight function is a ^-generalisation of a 
classical weight function. In [7] it is shown that the method of Borodin leads to (/-discrete 
Painleve equations. Can one obtain the q-discrete Painleve equations from an orthogonal 
polynomial approach? 
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